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Abstract 

The early E < 511 keV gamma-ray emission from novae depends critically on the ^*F(p,a)^^0 
reaction. Unfortunately the reaction rate of the ^*F(p,a)^^0 reaction is still largely uncertain due 
to the unknown strengths of low-lying proton resonances near the ^*F+p threshold which play 
an important role in the nova temperature regime. We report here our last results concerning 
the study of the d(^**F,p)^^F(a)^^N transfer reaction. We show in particular that these two low- 
lying resonances cannot be neglected. These results are then used to perform a careful study of 
the remaining uncertainties associated to the '^®F(p,q)'^^0 and '^®F(p,7)^®Ne reaction rates. 



Key words: Transfer reaction, Radioactive beam, DWBA, Spectroscopic factor, Nova nucleosynthesis. 
PACS: 26.60.Je, 21.10.Jx, 26.30.+k, 27.20+n 



Preprint submitted to Elsevier 



8 February 2008 



1. Introduction 




Gamma-ray emission from classical novae is dominated, during the first hours, by 
positron annihilation resulting from the beta decay of radioactive nuclei. The main con- 
tribution comes from the decay of ^^F (half-life of 110 min) and hence is directly related 
to ^^F formation during the outburst [l|, A good knowledge of the nuclear reaction 

rates of production and destruction of ^^F is required to calculate the amount of ^*F 
synthesized in novae and the resulting gamma-ray emission. The relevant reactions for 
the i^F production ["0(p,7)^SF and "0(p,a)"N] have been recently studied 0, S, i] 
and point to a lower production of ^^F. The relevant rates for the ^^F destruction are 
the i8F(p,7)i9Ne reaction {Q = 6.411 MeV) and the much faster i8F(p,a)i^0 reaction 
{Q — 2.822 MeV). In order to calculate these destruction rates, it is important to de- 
termine the properties of the ^^Ne states near the "'^^F -|- p threshold. Therefore, an 
extensive series of experiments like (p,p) elastic scattering direct (p,a) measure- 

ist and indirect transfer reactions to populate levels in ^^Ne and 
Il6l.ll7l| have been performed. However, its rate remains poorly known at nova 
temperatures (lower than 3.5x10^ K) due to the scarcity of spectroscopic information 
for levels near the proton threshold in the compound nucleus ^^Ne. This uncertainty is 
directly related to the unknown properties of the first three levels in ^^Ne: E^, — 
6.419 MeV, (3/2+); 6.437 MeV, (1/2") and 6.449 MeV, (3/2+) (see Fig 1) and following 
Utku et al. [14| we will assume the previous spin and parity assignments even though 
they have not been measured directly. The tails of the corresponding resonances at Er 
— 8, 26 and 38 keV, respectively, can dominate the -^*F(p,a)^^0 astrophysical S-factor 
in the relevant energy range 0]. As a consequence of these nuclear uncertainties, the ^*F 
production in nova and the early gamma-ray emission was estimated to be uncertain by 
a factor of r:!300 0. 

In order to estimate the level position and spectroscopic factors of the relevant states 
in ^^Ne, we performed shell-model calculations with the OXBASH code [l3|. It was 
not possible to reproduce experimental levels at energies about 6-7 MeV, motivating 
hence an experimental approach. Due to the very low Coulomb barrier penetrability for 
the relevant resonances in -'^^Ne, a direct measurement of their strength is at present 
impossible. Hence, we have used an indirect method, the d(^^F,p)^^F transfer reaction, 
aiming at the determination of the one-nucleon spectroscopic factor (S) in the analog 
levels of the ^^F mirror nucleus. Assuming the equality of the spectroscopic factors in 
analog levels, it is possible to calculate the ^^Ne proton widths through the relation 
Tp = S X Fg p 1 1 9l . |20| where Fg.p. is the single particle width. Preliminary results were 
presented in [15[, here we provide a detailed data analysis. 

This paper is organized as follows. The experimental method is detailed in Sections 2 
and 3. The data analysis is presented in Sections 4 and 5 and the remaining uncertainties 
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Fig. 1. Level scheme for and ^Ne above the ^F + p threshold. Dotted lines connect analog states 
from Il4ll . Dashed-dotted lines in ^^F represent the levels studied in the present work. 

concerning the ^*F(p,a)^^0 and ^*F(p,7)^^Ne reaction rates are studied in the concluding 
Section 6. 



2. Experimental method 
2.1. Experimental set-up 

The experiment was performed at the CYCLONE RIB facility at the Centre de Recher- 
ches du Cyclotron, UCL, Louvain-la-Neuve, Belgium. We used a 14-MeV "'^^F radioactive 
beam which was produced via the ^^0(p,n)^^F reaction, chemically extracted to form 
CHa^^F molecules, transferred to an ECR source, ionized to the 2+ state and then 
accelerated to the relevant energy [iH. The production of "'^^F was made in the batch 
mode and an average of 2.2x10^ ^*F ions per second on target was delivered for a total 
of 40 hours representing 15 ^^F batches. The beam contamination from its stable isobar 
was evaluated 20 hours after the delivery of the last batch from the counting rate on 
the target. The observed counting rate was in fact compatible with the ^^F radioactive 
decay and indicated a ^^0 / ^*F ratio less than 10"'^. 

The ^*F beam bombarded 100 /ig/cm^ deuterated polyethylene (CD2) targets. These 
targets were produced by polymerization and their thickness and homogeneity were 
checked by measuring the energy loss of a particles. The target stoechiometry and con- 
tamination were studied at the Orsay ARAMIS accelerator [22| after the d(^^F,p)^^F 
experiment. Elastically scattered protons were detected with a silicon detector at a labo- 
ratory angle of 150° for an incident proton energy of 2.8 MeV 2^, IJ]. It was found that 
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the stoechiometry was not modified during the experiment as expected from the low ^^F 
beam intensity. 

The products of the d(^^F,p)^^F reaction were detected using the multi-strip silicon 
detector LEDA [2^. Two configurations named LEDA and LAMP were used. These 
detectors are composed of 8 and 6 sectors, respectively, each divided into 16 radial strips. 
The energy calibration of the 224 strips was performed with a calibrated 3a-source 
(^^^Pu, ^''^Am and ^^"'Cm) whereas the time-of-flight calibration was performed with a 
precision pulser. 




LAMP 

LEDA 

Fig. 2. The experimental set-up is shown with ^^F ions impinging on a CD2 target. Protons were detected 
in LAMP while ^^N or a-particles from ^''F decay were detected in coincidence in LEDA. See Fig 1 
in [l5| for a perspective view of the set-up. 



The experimental set-up is shown in Fig 2. The protons from the d(^^F,p)^^F reaction 
were detected in LAMP positioned 8.7-cm upstream from the target. It corresponds to 
laboratory angles between 110° and 157°, i.e. center-of-mass angles between 12° and 
44° for the two 3/2+ levels of astrophysical interest. Since these two levels in -'^^F are 
situated above the a-threshold at 4.013 MeV, they decay through ^^F* a+^^N. The 
decay products were detected in LEDA positioned 40-cm downstream from the target, 
which corresponds to laboratory angles between 7° and 18°. The positions of the detectors 
were determined by Monte-Carlo simulations: for LAMP, the angular range and detection 
efficiency were maximized and for LEDA, the coincidence efficiency between protons 
detected in LAMP and ^^N detected in LEDA was optimized. 

Rutherford elastic scattering of ^^F on the carbon contained in the target, detected in 
LEDA, provides the (target thickness) x (beam intensity) normalization. We used the 
strip at 15.7° in LEDA for which the solid angle is maximum and less sensitive to the 
exact position of the beam on the target. Furthermore, for this strip the elastic scattering 
peaks of ^*F and ^^C are well separated. The uncertainty on the normalization is mainly 
due to the position of the beam (± 2 mm) and is estimated to be of 7%. 

Data were collected in event-by-event mode where the multiplicity, the angle, the 
deposited energy and the time of flight relative to the cyclotron radio-frequency were 
recorded, allowing an off-line analysis of single and coincidence events. 
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2.2. Data reduction 

Thanks to the low ^^F energy {Ec.m. = 1-4 MeV), below the Coulomb barrier of the 
+ ^^C and ^^F + d systems, only a few channels are opened: d(^^F,a)^^0, d(^^F,p)i9F, 
d(-'^*F,pa)-'^'''N, d(^^F.n)^^Nc, and d(-'^*F,na)^'''0. Single and coincidence events have been 
analyzed but the coincidence condition between LAMP and LEDA (see Figure 3) al- 
lows to distinguish easily the events from these reactions. The two-body properties of 
the d("'^*F,Q;)^^0 reaction can be observed as a linear correlation between the detected 
energies in LAMP and LEDA, which is not the case for the three-body d(^*F,pa)^^N 
reaction. We performed a further check of the reaction identification by verifying that 
the reaction products were (not) detected in the same reaction plane, indicating a two 
(three)-body reaction. 
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Fig. 3. Ei^AMP X l^LEDA coincidence spectrum showing several regions corresponding to the d(^**F,p)^''F 
(solid lines regions) and d('^*^F,a)^''0 (dashed lines regions) reactions. Label 1 indicates the d(^*F,pa)'^^N 
reaction while label 2 indicates the d(^**F,p)^''F reaction for ^^F excited states below the a— threshold. 
Label 3, 4 and 5 correspond to different ^^O excited states for the d(^*F,a)i^O reaction. 



Kinematical bands (energy versus strip number in LAMP) corresponding to the events 
in Fig 3 are plotted in Fig 4 (right). For comparison, kinematical bands from single events 
are also plotted in Fig 4 (left). The same patterns as for coincidence events are observed 
but, in addition, since there is no coincidence condition, some kinematical bands are 
completed (region 2) or new (between 2 and 4 MeV). 

Protons from the d(^^F,pQ;)^^N reaction are detected in LAMP in coincidence with 
a-partick; or ^''N in LEDA. Both type of coincidences can be separated using the time- 
of-flight information in LEDA, but since the p-a coincidences represented only 7% of the 
total coincidence events, only the p-^^N coincidences were analysed. After selection of the 
p--'^^N coincidences, the excitation energy of the decaying ^^F levels can be kinematically 
reconstructed from the energies and angles of the detected protons and the known beam 
energy. From the ^^F excitation energy spectra, and after identifying the populated levels, 
the differential cross section in the center-of-mass system was calculated following the 
formula: 
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Fig. 4. Kinematical bands in LAMP for single and coincidence events. Angles are given in arbitrary 
units. Labels have the same meaning as in Fig 3. The different slopes for the kinematical bands (region 
4 and 1) is the signature of different reactions [d(l^F,a)l^O and d(^*F,p)'^®F, respectively]. 

\dCl J ^ ' NisF Nd AO 

where Np is the number of detected protons in a given LAMP strip corresponding to the 

angle 9, Nis p and Nd are the number of incident -'^^F ions and the deuteron content of 
the target per unit area, respectively, and AfJ is the geometrical solid angle. In the case 
of coincidences, Ai7 is multiplied by the p-^^N coincidence efficiency vi^hich is determined 
by Monte Carlo simulations considering an isotropic angular distribution in the a-^^N 
center-of-mass system for the a-decay of the ^^F. However for the two 3/2+ levels, the 
a-particle is emitted with an £ = 1 angular momentum. Monte-Carlo simulations using 
a-particle angular distribution corresponding to different sub-magnetic population, gave 
an uncertainty related to the normalisation of about 15%. 



3. Spectra and peak deconvolution 

3.1. i^F Spectra 

The reconstructed ^''F excitation energy spectra are shown in Figure 5 for coincidence 
(bottom) and single (top) events after applying the alpha calibration and before the 
internal calibration (see the end of the subsection for more details). Owing to the high 
coincidence efficiency (e = 70%), the coincidence and singles counting statistics were 
comparable. Both spectra correspond to all the events collected in the 6 sectors of the 
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LAMP detector without any correction for the beam position on target. Because LAMP 
is very close to the target, this reconstruction is very sensitive to the beam position on 
the target, which is known within ± 2 mm. Hence the angular position of the strips 
on a same ring varies, inducing a different reconstructed excitation energy for the same 
detected proton energy. This, in turn, leads to a deterioration of the resolution when 
the contribution of the 6 LAMP sectors are directly summed without correction. To 
compensate this effect, the 6.5 MeV peak obtained in coincidence for each LAMP sector 
is fitted and each spectrum is shifted relative to the mean position of the 6 sectors. 

For the coincidence events, the vertical lines in Figure 5 (bottom) present the position 
of the ^^F levels [l^ populated with low transferred angular momentum {i < 2). Even if 
the energy resolution (FWHM « 100 keV) is not sufficient to separate the 20 ^^F levels 
that are represented, the two 3/2+ levels of interest in ^^F at = 6.497 and 6.528 MeV 
(the analogs of the 3/2+ levels in ^^Ne) are well isolated from the other groups of levels. 
Unlike for coincidence events, single events (Figure 5 top) are directly selected from a 
time-of-flight versus energy spectrum in LAMP, which does not allow to separate proton 
events from a-particle events coming from the d(^*F,a)^^0 reaction. For ^®F excitation 
energies higher than about 5 MeV, one can observe the same structure for the single 
and coincidence spectrum. However, it should be noted that the 7.2-7.3 MeV peak is less 
populated than in the coincidence spectrum because a higher energy threshold for the 
selection of single events is used to reduce electronic noise. For lower excitation energy 
(below the a-threshold), ^^F levels are observed down to the ground state and the vertical 
lines represent all the existing "'^^F levels. 

We performed an internal excitation energy calibration using the peaks corresponding 
to the well isolated ^^F levels at = 2.780 MeV, 5.106 MeV as well as the doublet E^ 
— 7.262 + 7.363 MeV. In case of the unresolved peak at 7.3 MeV, its doublet nature 
was checked by performing a fit with two components where the position of each compo- 
nent was left as a free parameter. The experimental energy difference obtained {AE^ = 
107 keV) is in very good agreement with the energy difference from the literatu re ( AEx 
= 101 keV) and the hypothetical presence of a third state in this energy region [l7| was 
not needed to describe the data. 



3.2. Peak deconvolution 

The "'^^F excitation energy spectrum obtained after the internal calibration is shown in 
Figure 6 as well as a fit to the data between E^ = 6150 and 6950 keV. The fit includes the 
contribution of two levels for the 6.5 MeV peak {E^ = 6497-1- 6527 keV), three levels for 
the 6.25 MeV peak {E^ = 6255 + 6282 -I- 6330 keV) and three levels for the 6.9 MeV peak 
{Ex — 6787 -f 6838-1-6891 keV). The free parameters of the fit are the amplitude of each 
gaussian and a width common to all the levels. The level energies are fixed. The excess of 
counts below the 6.5 MeV peak is interpreted as the broad l/2~ level E^ = 6536 keV 
(r = 245 keV) for which only the intensity is left as a free parameter. 

For this broad ^^F 1/2^ level, we used the values of the energy and width derived 
recently [27j from a reanalysis of previous data [2^. Compared to previous adopted 
values, the width is similar (245 keV instead of 280 keV) but the energy is about 100 keV 
higher {E^ = 6536 keV instead of E^ = 6429 keV). It should be noted that the fit obtained 
with the new parameters is better (x^/ndf = 1.19) than with the previous parameters 
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Fig. 5. Top: Reconstructed ^^F spectrum for single events after the alpha calibration. Levels down to the 
ground state are populated. Contribution from the d(^*F,«)'^^0 reaction is shown (see text for details) 
Bottom: Same for coincidence events, but covering only excitation energies above the a-particle emission 
threshold. The doublet of 3/2"*" astrophysical levels around 6.5 MeV is well isolated from the other levels. 
Vertical lines show ^^F levels (see text). 



(x^/ndf — 1.54, see Figure 1 in \29\). Furthermore, a fit with the same conditions but 
with the energy of the broad l/2~ level as free parameter naturally gives an energy of 
Ex = 6529 ± 13 keV, compatible with the new value from [27] . 

The solid line on Figure 6 shows the best fit and the dashed lines represent the contri- 
bution of each level separately. For the 6.5 MeV peak composed of the two 3/2+ levels, 
the high-energy level is favored with a relative contribution of 5%-95% for the = 6497 
and 6527 keV level respectively. A recent experiment using the same reaction but a higher 
bombarding energy found an opposite proportion for the two 3/2+ levels [17| . In both 
cases, an internal energy calibration of the ^^F excitation energy was used. The differ- 
ences are that our energy calibration is an interpolation based on experimental peaks 



;xpe i 

surrounding the 3/2+ levels of astrophysical interest whereas in Ref. [l7| the calibration 
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Fig. 6. Reconstructed F coincidence spectrum after applying the internal energy calibration. The solid 
line represents the best fit of the spectrum and the vertical lines indicate the involved levels. Each level 
contribution is shown in dashed lines. See text for details. 

is an extrapolation done for several laboratory angles, from peaks at lower energies than 
the levels of astrophysical interest. 

However, since the two levels are not experimentally resolved, different relative con- 
tributions can give similar ^ and hence the nominal relative contribution derived above 
should be taken with care. Hence we performed a statistical study aiming at the deter- 
mination of the upper limit of the Ex — 6497 keV level contribution. A careful study 
of the systematic errors (target thickness and detector dead layer, beam energy, energy 
loss in the target and in the detector dead layer, beam centering and detector position) 
associated to the reconstruction of the excitation energy, showed that the most important 
uncertainty is related to the precise position of the LAMP detector. Hence three uncer- 
tainties were taken into account: the error associated to the internal energy calibration 
(k, 3 keV for ~ 6.5 MeV), the error associated to the precise position of LAMP after 
the internal energy calibration (« 5 keV for ±1 mm) and the deconvolution procedure 
of the 6.5 MeV peak. The l-tr upper limit for the contribution of the E^. — 6497 keV 
level for each effect taken separately is 41%, 34% and 44% respectively. When summed 
quadratically, we obtain a contribution for the E^ — 6497 keV level of 5^5^%- It should 
be noted that there is a 7/2 level of unknown parity at Ex = 6554 keV lying just above 
the 3/2+ doublet. We estimated its contribution to be less than 5%. 

From the previous analysis, one cannot exclude that the two 3/2"*" levels are populated 
in the same proportion. However one expects that these two levels have different structure 
properties as it is suggested by an electron inelastic scattering measurement on ^^F [sot . 
Such an experiment gives access to the tranversal and longitudinal form factors [i^^(Q')] 
as a function of the transfered momentum . For such an inelastic scattering, the form 
factor depends on the initial and final nuclear states and hence of the nuclear structure 
of the observed levels. The longitudinal form factor of the two 3/2+ levels are different by 
more than one order of magnitude, indicating a very different nuclear structure. It is then 
not surprising that one of the two levels would be preferentially populated through the 
d("'^^F,p)"'^^F transfer reaction, since transfer reactions are very sensitive to the nuclear 
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structure. 



4. DWBA and compound nucleus components 

4.1. DWBA parameters 

DWBA calculations were performed with the FRESCO code ^l[ in Finite-Range. For 
testing the dependency of the results versus the DWBA parameters, two sets of optical 
parameters were used: A+A' and B+B' (see Table 1). The optical potentials are written 
as usual: 



Vopt Vb - Vo/(r, To, ao) 



Wf{r, rw,aw) - 4Wd — /(r, rD,aD) 
ar 

, Vs.o 



~] ^■^-rfi'^^ rs.o.,as.o.) 
r \mTrC / ar 

where f(r,ri,ai) is a potential well of Wood-Saxon shape with r, and a^ being the radius 
and the diffuseness of the interaction potentials. Vc, Vq, W, Wd and Vg.o are the 
Coulomb, the volume, the volume absorption, the surface absorption and the spin-orbit 
well depths, respectively. Usually the optical parameters are determined by fitting elastic 
scattering data of the entrance and exit channels. Such data do not exist for the d+^*F 
channel, so here the parameters correspond to other nuclear systems of similar masses 
and similar energies. The first set of parameters A+A' comes from a compilation [s^ ]. 
whereas the second set of parameters B+B' comes from a similar neutron transfer reaction 
"'^^F(d,p)^°F which was studied at the same center-of-mass energy (s^- In that study, 
11 angular distributions of ^'^F excited states are well described by the optical potential 
B+B'. For the Finite-Range calculation, the potential describing the deuteron is of Reid 
Soft Core type [34|. The depth of the volume potential well for the transferred neutron 
is automatically adjusted to reproduce its separation energy in ^^F. 

Table 1 

Optical potential parameters used for the DWBA analysis. Potentials A+A' are from fs^] and B+B' are 
from [331 . Notations are explained in the text. 



Optical 
potentials 



A 
B 



A' 
B' 



rc Vo ro ao W rw Wd to ajj N s.o. rs.o. as.o. 

(fm) (MeV) (fm) (fm) (MeV) (fm) (fm) (MeV) (fm) (fm) (McV) (fm) (fm) 



+ d entrance channel 
1.3 80.1 1.1 0.972 — — — 14.8 1.6 0.652 
1.3 60.0 1.5 0.6 20 1.5 0.6 — — — 



1.3 47.45 1.185 0.721 
1.3 50.0 1.3 0.5 



p exit channel 

— 7.5 0.942 0.568 5.1 1.042 0.488 

— 8 1.3 0.5 — — — 



It should be noted that even if the energy available in the center-of-mass system 
(1.4 MeV) of the d(^^F,pQ;)^^N reaction is lower that the Coulomb barrier (Ec = 
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2.6 MeV), the emitted protons from the two 3/2+ levels of astrophysical interest are 
slightly above the Coulomb barrier of the exit channel {Ep = 3.1 MeV and Ec = 

2.7 MeV). Hence the term of "sub-Coulomb transfer" does not seem to be justified 
for the conditions of this reaction. We have checked this by computing a DWBA calcula- 
tion in which the nuclear potentials were switched off, only leaving a Coulomb potential. 
This calculation was unable to reproduce the shape of the angular distribution for the 
two 3/2+ levels. 



4.2. Compound nucleus component 



Due to the restricted angular range {0c.m. < 50°), we used the 9/2+ E^ = 2780 keV 
"'^^F level observed in our data to estimate the compound nucleus component. Its spin and 
parity make it difficult to be populated by a transfer reaction due to the high transferred 
angular momentum £ = 4 which corresponds to a 9.2 MeV centrifugal barrier, to be 
added to the Coulomb barrier. The DWBA calculation for this level with a maximum 
spectroscopic factor of 1 (see Figure 7) fails to describe the data. Two-step processes such 
as the excitation of the 3+ state at 937 keV in ^^F followed by the transfer of a neutron 
into an £ = 2 orbital were also computed with the coupled channels code CCZR [3l|. The 
contribution of these two-step processes was at most a few /ib/sr and can not explain the 
angular distribution of the 9/2+ level. The Hauser-Feshbach contribution was calculated 
with the HSFB [3^ code taking into account the open channels from the compound 
nucleus ^'^Ne decay at our bombarding energy and normalized to the E^ = 2780 keV data. 
Applying this normalization to the 5/2+ E^ ~ 5106 keV ^^F level describes successfuly 
the data which does not seem to be described by a dominant direct reaction mechanism 
(see Figure 7). 



5. Discussion 



The angular distributions of single data for the E^i^^F) — 2780 and 5106 keV levels and 
of coincidence data for the doublet E^(^^F) = 6497 and 6527 keV and the two E^(^^F) = 
7262 and 7364 keV levels are shown in Figure 7. While the angular distributions for the 
two levels at E^ — 2780 and 5106 keV are quite flat, supporting a compound nucleus 
type reaction mechanism, the angular distribution of the 3/2+ doublet shows a strong 
variation of about one order of magnitude in a small center-of-mass angular range, which 
is typical of a direct reaction mechanism. The spectroscopic factors obtained from our 
analysis are displayed in Table 2. 



5.1. Levels at E^ = 6.497 and 6.527 MeV m ^^F 



The two 3/2+ levels can be populated by both I — Q and 1 = 2 transfers, however the 
fit of the DWBA calculations to the data indicates that only the t = component is 
contributing (see Figure 7) at such low center of mass energy. At much higher incident 
energy, a £ = 2 component was requested to reproduce transfer data 17|. Taking into 
account the compound nucleus contribution as explained above, we obtain for ^ = a 
spectroscopic factor of S'(6.497 -I- 6.527) = 0.17 (to be compared with S = 0.12 from a 
recent experiment [l3l)- This value supersedes our previous value of 5(6.497 -I- 6.527) = 
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Fig. 7. Angular distributions for the i^F levels at 2.780, 5.106 MeV, the doublet at 6.497 and 6.537 MeV, 
and the levels at 7.262 and 7.364 MeV, populated by the d(^^F,p)^^F transfer reaction. Dashed lines 
represent the compound nucleus component and solid line are the results of the direct contribution 
calculated using a Finite-Range DWBA analysis and the B+B' optical potentials. 



0.21 [l5|: the difference comes from the new Finite-Range DWBA calculations and the 
inclusion of the Hauser-Feshbach contribution. The fact that the two sets of optical 
parameters A+A' and B+B' give very similar angular distributions and spectroscopic 
factors makes us to feel confident about the present results. The analysis of the single data 
events for this 3/2+ doublet allowed us to derive the corresponding angular distribution 
and to compare it with the coincidence one. They agree well, both in the shape and in the 
normalisation, hence giving strong confidence in the coincidence solid angle determined 
by Monte-Carlo simulations. 

The uncertainties associated to the spectroscopic factors have been evaluated to be 
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Table 2 

Spectroscopic factors for the ^^F levels analyzed in DWBA. 



E^(19F) [MeV] 


J" In 


C^S 


6.497 1 








3/2+ 


0.17 


6.527 J 






6.536 


1/2- 1 


< 0.33 


7.262 


3/2+ 


0.15 


7.363 


1/2+ 


0.22 



of 25%. The major source of uncertainty comes from the DWBA method (20%) and is 
related to the choice of the optical parameters in the entrance channel which cannot 
be constrained by experimental data. Another source of uncertainty (15%), which has 
already been discussed, arises from the angular distribution of the emitted a-particle in 
the ^^F center-of-mass frame. The uncertainty on the normalization contributes to 7% 
(Section 3). 

5.2. Level at ^ 6.536 MeV in ^^F 



Due to its large total width F ~ 245 keV [27[, this 1/2 level can play an important role 



in the ^*F(p,a)^^0 reaction since it covers the whole Gamow peak at novae temperatures. 
But its large width also makes it difficult to be observed. However, making the hypothesis 
that all the excess of counts below the 6.5 MeV peak comes from the 1 /2~ level, an upper 
limit to its contribution could be derived. This hypothesis does not seem unrealistic 
because no counts are observed in the coincidence spectrum in the = 5.8 MeV region 
(see Figure 5 bottom), where no ^^F levels exist. Hence this seems to indicate that the 
excess of counts below the 6.5 MeV peak has a physical origin. 

The number of counts corresponding to the population of the 1/2^ level is obtained 
from the global fit presented in Figure 6 where this level is described by a lorentzian shape 
whose amplitude is a free parameter. Then, assuming a £ = 1 transfer for populating this 
level, one obtains an upper limit for the spectroscopic factor of S'(6.536) < 0.33. The 
difference with respect to our previous value (£'(6.536) < 0.15 [15]) comes from the new 
position and width of this level [27j] as well as from the DWBA calculation, now made in 
Finite-Range. 



5.3. Levels at E^ = 7.262 and 7.364 MeV in ^^F 



The experimentally unr esolved ^^F levels at 7.262 and 7.364 MeV have a 3/2+ and 
1/2+ spin and parity, respectively, and both can be populated bya£=:0or£ = 2 
neutron transfer. However the DWBA calculations fitted to the data showed that only 
the £ — is contributing (see Figure 7). These two levels have the highest excitation 
energy which can be achieved with the beam energy used here, and correspond to the 
less energetic protons detected {Ep < 1 MeV). Part of the kinematical band of these 
two levels, which are not experimentally resolved, is cut by the electronic threshold. 
Hence only the five outer strips of LAMP were considered in the analysis. The extracted 
angular distribution does not allow to deduce the relative contribution of the two levels 
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since the differential cross section has the same shape due to the same £ = transfer 
but a different normalisation (a factor of 2) due to the spin factor. However a two-peak 
fitting analysis of the 7.3 MeV excitation energy peak indicates a relative contribution 
of 57%-43% for the 7.262 and 7.364 MeV levels, allowing to obtain a set of spectroscopic 
factors S'(7.262) = 0.15 and S'(7.364) = 0.22. The compound nucleus component which 
has been taken into account has a negligible cross section (« 0.1 mb/sr). 

5.4. Analog level of the E^{^^Ne) = 7.076 MeV 

Up to now, the analog level in ^^F of the ^^F(p,a)^^0 resonance at Er — 665 keV 
[E^i'^^'Ne) = 7.076 MeV, J'' = 3/2+] is still not identified. A gamma-ray measure- 
ment via the ^^N(a,7)i^F reaction found a level at 7.101 MeV in ^^F with a width 
Fq, = 28 keV. The spin and parity was deduced to be 3/2+ by observing gamma-ray 
transitions to a 3/2~ level. Assuming equality of reduced a-widths between analog states. 
Butt et al. [33 find r„ « 30 keV for the ^^Ne level at E^ = 7.076 MeV, which is in good 
agreement with the value of Fq, = 24 keV from [Ot] . 

However, the assignation of a 3/2+ spin and parity to the E^C^F) — 7.101 MeV level 
has been questioned by Fortune and Sherr [38*1 who pointed out that this level could be 
a doublet of spins 1/2-5/2 or 1/2-7/2. Furthermore, these authors predict the position 
of the analog level at higher energy around E^ — 7.4 ±0.1 MeV. 

Assuming here that the 7.101 MeV level in ^^F is the analog of the E^C^^e) = 
7.076 MeV and assuming equality between dimensionless reduced widths [8^(^^F) = 
6p(^^Ne)], we obtain 0^(^^F) — 0.14. This value is comparable to the one from the two 
3/2+ levels of astrophysical interest (8^j — 0.1) and should hence appear in our spectra 
as a peak of similar intensity as the 6.5 MeV peak, but our data (Figure 5, bottom) 
do not show any evidence of a peak around 7.1 MeV. Furthermore it has been shown 
recently that the spin of this resonance was unlikely to be 3/2+ (27j . 

The 7.3 MeV group in ^^F is composed of two states at 7.262 MeV (3/2+) and 
7.368 MeV (1/2+) whose spectroscopic factors are determined here to be = 0.15 and 
0.22, respectively (see Table 2). The spectroscopic factor value for the 3/2+ state is very 
close to the one of the level at 7.076 MeV in ^^Ne, and seems then to indicate that the 
E^C^F) = 7.262 MeV is a good candidate for the analog level of the 7.076 MeV in ^^Ne. 
However, it should be noted that the total width of the 7.262 MeV level in -'^^F calculated 
from the Ea;C^'Ne) = 7.076 MeV is F = Fq, « 28 keV which seems incompatible with the 
measured value F < 6 keV [2^. 

5.5. Proton widths in ^^Ne 

The obtained neutron spectroscopic factors are for ^^F levels whereas the proton widths 
of astrophysical interest correspond to their analog levels in ^^Ne. Following Utku et 
al. Hi], we assume that the two 3/2+ levels E^ = 6497 and 6527 keV and the 1/2" level 
E^ = 6536 keV populated in ^^F are the analog states of the 3/2+ levels E^(^^l<le) = 6419 
and 6449 keV, and of the Ex{^^Ne) = 6437 keV, respectively, even if the spins and parity 
of ^^Ne states has not been determined experimentaly. Unlike other analog levels in 
i9p^i9]\jg where the assignation comes from the fact that the levels are populated with 
comparable statistics by mirror reactions, the case of the two 3/2+ and l/2~ levels only 
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relies on the similarity of their energy position and of their total width, respectively. 
The two 3/2+ levels are both separated by 30 keV in ^^F and ^^Ne at similar excitation 
energies, whereas for the 1/2" level r„(i9Ne) = 220 keV « TaC-^F) = 245 keV. 

It is of common use in nuclear astrophysics to assume equality of spectroscopic factors 
between analog states. Hence for the 3/2+ states in ^^Ne, we obtain a proton spectroscopic 
factor S = 0.17 and for the l/2~ state, we obtain an upper limit S < 0.33. However, 
the uncertainty associated to this practice is estimated to be a factor of two from the 
comparison of the values of neutron and proton spectroscopic factors for all the analog 
states of the sd shell nuclei (Figure 8). 
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Fig. 8. Comparison of neutron and proton spectroscopic factors for analog levels of the sd shell nuclei. 
Data comes from [l9| . 



The proton widths are calculated using the following relations: 

Tp^C^STs.p. with r,.p. = — |7^,.p.(s)|'P^, 

where Ts.p. is the single-particle width, s is the interaction radius (s = 4.5 fm) and 
TZs.p. is the proton radial wave function in the p+^*F system calculated with classical 
parameters for the Wood-Saxons nuclear potential: (r, a) — (1.25, 0.65) fm, where r and 
a are the radius and diffuseness of the potential. As compared to direct determination 
of the proton width, a statistical study of sd shell nuclei showed an agreement between 
direct and indirect method of about 30% [s^ . Using the above prescription we obtain for 
the 3/2+ doublet in ^^Ne, rp(6.419) = 3.4 x lO'^^ keV and rp(6.449) = 1.9 x lO'^^ keV 
assuming that the total strength of the doublet is put on each level respectively. As far as 
the levels above = 1 MeV are concerned and if we assume that the levels Ex{^'^Y) = 
7.262 MeV and E^{^'^'^e) = 7.238 MeV are analogs [13 we obtain rp(7.238) = 10.3 keV 
which is larger than the adopted upper limit (Fp < 4 keV [31 )■ For the analog state of 
the E^{^^¥) = 7.364 MeV 1/2+ level, no 1/2+ levels in ^^Ne have been observed so far 
at an excitation energy around 7.3 MeV. 
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6. CONCLUSIONS: REMAINING UNCERTAINTIES 



Our analysis showed that the two expected low-lying 8 keV and 38 keV proton reso- 
nances {Ex = 6419 and 6449 keV in ^^Ne) must be considered for the calculation of the 
^*F(p,a)"'^^0 reaction. Here we focus on the remaining uncertainties before firm reaction 
rates can be computed for the ^*F(p,a)^^0 and ^*F(p,7)^^Ne reactions. 

Due to the small energy difference (30 keV) between the two 3/2+ levels in ^^F {E^ — 
6497 and 6527 keV) and the two resonances in ^^Ne {Er = 8 and 38 keV), the assign- 
ment of analog pairs is not straigthforward and one cannot discard the possibility of an 
inversion. Hence the identification of the main contributing ^^F level is not so crucial 
for the determination of the associated proton widths in ^^Ne. So two cases will have to 
be considered, the whole strength at the Er = 8 keV resonance or at the E,. = 38 keV 
resonance. 

The second uncertainty arise from interference effects between the two Er = 8 and 
38 keV resonances and the 3/2+ resonance at Er = 665 keV (corresponding to the ^^Ne 
level Ex = 7.076 MeV). The sign of these interferences is totally unknown and for the 
destructive case, the S-factor would become extremely small between 100 and 300 keV, 
which corresponds to Gamow peak energies for temperatures achieved in novae (50 < 
Tg < 350). As usual, the notation = y corresponds to the temperature y x 10"^ K. 

Another uncertainty on the reaction rates comes from the poor knowledge of the ^^Ne 
alpha widths, Fq, for the Er = 8 and 38 keV resonances. Since there are no measurements 
of these F^ in ^^Ne, they are usually deduced from the analog states in -^^F. However 
differences as large as a factor 10 have been observed for alpha widths of analog levels 
for the ^^F / ^^Ne mirror nuclei Calculations of constructive and destructive inter- 
ferences for different values of alpha widths showed a strong impact on the ^*F(p,Q;)^^0 
reaction rate 41 1. 

Finally the energy of the two Er = 8 and 38 keV resonances in ^^Ne is not known 
better than 6 keV [lj|. This effect is far less important than the one of the unknown 
alpha widths because the position of the resonances are below the energies corresponding 
to the Gamow peak for temperatures relevant to ^^F production in novae. 

In conclusion we reported here on our latest results concerning the one-nucleon transfer 
reaction d(^^F,pQ;)"'^^N used to study -'^^F analog states of the astrophysically important 
states in ^^Ne. A DWBA analysis of the data showed that the two 3/2+ low-energy 
resonances in ^^Ne, Er ^ 8 and 38 keV have a relatively large spectroscopic factor (5 
= 0.17) and hence cannot be neglected in calculating the ^*F(p,a)-'^^0 and ^*F(p,7)-'^^Ne 
reaction rates. We made a careful analysis of the impact on the reaction rates of remaining 
uncertainties. This includes interferences between the low-lying states and the 3/2+ £",. = 
665 keV resonance, the unknown a-widths in ^^Ne and the possibility of inversion for 
the analog levels between ^^F and ^^Ne. A new challenging experiment to investigate 
the ^'*F(p,a)^^0 reaction at energies below 0.6 MeV in the center-of-mass will be soon 
performed at the Louvain-la-Neuve CYCLONE facility. 
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